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ABSTRACT 

We present results from a deep Suzaku observation of the black hole in LMC X-l, 
supplemented by coincident monitoring with RXTE. We identify broad relativistic re- 
flection features in a soft disc-dominated spectrum. A strong and variable power-law 
component of emission is present which we use to demonstrate that enhanced Comp- 
tonisation strengthens disc reflection. We constrain the spin parameter of the black 
hole by modelling LMC X-l's broad reflection features. For our primary and most 
comprehensive spectral model, we obtain a high value for the spin: a* = 0.97^g^3 (68 
per cent confidence). However, by additionally considering two alternate models as a 
measure of our systematic uncertainty, we obtain a broader constraint: a* = 0.9Tlto 25- 
Both of these spin values are entirely consistent with a previous estimate of spin ob- 
tained using the continuum-fitting method. At 99% confidence, the reflection features 
require a* > 0.2. In addition to modelling the relativistically broadened reflection, we 
also model a sharp and prominent reflection component that provides strong evidence 
for substantial reprocessing in the wind of the massive companion. We infer that this 
wind sustains the ionisation cone surrounding the binary system; this hypothesis nat- 
urally produces appropriate and consistent mass, time, and length scales for the cone 
structure. 

Key words: accretion, accretion discs — black hole physics — stars: individual 
(LMC X-l) — X-rays: binaries 



1 INTRODUCTION 

LMC X-l was the first extragalactic black hole (BH) binary 
to be discovered; Cygnus X-l is the only other such persis- 
tent X-ray source with an O-giant companion that is located 
locally, i.e., within the Galaxy and Magellanic Clouds. LMC 
X-l is quite unusual in that it consistently maintains a sta - 
ble luminosity of -Lx/^Edd ~ 16 per cent l|Gou et all 120091 ') 
despite showing strong fluctuations in the rms amplitude of 
its power spectrum. 

Recently, the spirfl of LMC X-l was measured by mod - 
elling the thermal accretion disc emission (iGou et alj|200gft 
via the X-ray continuum-fitting technique (e.g., Zhang et al.1 
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1 Defined as the dimensionless parameter a* = cJ/GM 2 with 
a* | ^ 1, where M and J are the BH's mass and angular momen- 
tum. 



1997ft. Usi ng a primary sample of 18 RXTE PCA spectra, 



Gou et all estimated the spin to be a, — 0.921q'q7 



The principal alternative to continuum fitting is the re- 
flection (or Fe-line) method. Here, the breadth of relativisti- 
cally broadened reflection features generated in the accre- 
tion disc, most notably the prominent Fe Ka line com- 
plex, is used to determine spin (e.g.. iFabian et all 1 19891 : 
iBrenneman fe Revnold"sll2006l ; lMillerll2007ft . One advantage 
of the Fe-line method over X-ray continuum fitting is that it 
can be readily applied to measure the spins of stellar-mass 
and supermassive BHs alike. In the case of stellar-mass BHs, 
the Fe-line method has the further virtue that it is indepen- 
dent of BH mass and distance; however, its primary draw- 
back is t hat it relies upon a much fainter signal. Recent stud- 
ies (e.g.. iMiller et al.|[2009l : ISteiner et alj|201lft have made 
headway in achieving measurements via both techniques. 

Both methods make one fundamental assumption: that 
the accretion disc is truncated at the innermost stable cir- 
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cular orbit (ISCO). The radius of the ISCO is uniquely 
defined by the BH's mass and spin, growing with increas- 
ing mass, and shrinking with increasing spin. The ISCO- 
truncation assumption is presently at the forefront of the- 
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but see 



le et al.l 12009 . and references therein). At the 
same time, several empirical studies of BH binaries in disc- 
dominated states provide strong support for the existence of 
a stable inner-disc radius, including a compre hensive study 
analy sing decades of spectra of LMC X-3 jSteiner et al.l 
I2010T I and other studies of BH spectral evolu tion during out- 
burst (e.g.. iDone et al ] l2007l : lGilfanovll2010l ). 

Fortuitously, BH X-ray binary sources located in the 
LMC are in Suzaku's 'Goldilocks zone' for Fe-line observa- 
tions. By virtue of being located at ~50 kpc distance, even 
luminous BHs are faint enough that pileup effects are mod- 
es10, while at the same time these sources are bright enough 
to provide the signal-to-noise required for reflection spin 
measurements in a reasonable observation time of ~100 ks. 
By comparison, pileup in standard Suzaku observing modes 
is severe for outbursting Galactic sources, while Fe-line spin 
measurements for practical observing times would be photon 
starved for stellar-mass BHs beyond the Magellanic Clouds. 

LMC X-l is propitious for study because, despite being 
in an almost persistently disc-dominated state, its spectrum 
contains a relatively strong power-law component (that is re- 
quired to produce reflection). Tentative evidence for broad 
Fe Ka e mission (with equival ent width ~ 200 eV) was re- 
ported bv lNowak et al.1 (|200ll ). We unambiguously confirm 
the presence of a broad Fe line in LMC X-l, and we use 
it along with other reflection features to constrain the BH's 
spin. We thereby demonstrate that Fe fluorescence emission 
is present in the soft states of LMC X-l and we furthermore 
show that the strength of this feature is directly related to 
the strength of the Compton power-law component. 

Our Suzaku and RXTE observations of LMC X-l are 
described in Section [2] and the spectral models we employ 
are described in Section [3] Our analysis is presented in Sec- 
tion [4] followed by a discussion in Section[5l and we offer our 
conclusions in Section [6] 



2 OBSERVATIONS 

We observed LMC X-l for 130 ks using Suzaku from 2009 
July 21 through 2009 July 24 with the XIS detectors oper- 
ating in quarter-window mode. The data were reduced fol- 
lowing the XIS and PIN pipeline procedures. For the XIS, 
we found it necessary to increase the sisclean threshold to 
prevent contamination of the image corj_]. During the ob- 
servation, the source count rate was stable, varying from its 
average intensity by less than 5 per cent. Background spec- 
tra were obtained from observations of the Lockman Hole. 
In this and other reduction steps, we have been guided by 
the procedures described in iKubota et al.1 (|2010T l . 



2 iMiller et al. (2010) discuss pileup and its impact on line fitting. 

3 See http://www-utheal.phys. s. u-tokyo.ac.jp/ 



Suzaku's atti tude calibration w as improved using the 
aeattcor routincQ ( Nowak et alJl-ioill ). The innermost region 
of the point-spread function (PSF) suffered from moderate 
(~ 10 per cent) photon pileup. To ameliorate this problem, 
we used the utility pile-estimat^\ as a guide, and excised the 
innermost 30". Doing so, we retained « 70 per cent of the 
flux, while keeping the net pileup below 3 per cent. 

Three XIS units were used to collect the data: the two 
front-illuminated detectors, XIS-0 and XIS-3, and the back- 
illuminated detector, XIS-1. Spectra were binned to approx- 
imately half Suzaku's energy resolution and analysed over 
0.8-10 keV (XIS-0,3); 0.8-8 keV (XIS-1). Because of cali- 
bration defects, we omitted channels between 1.5-2.5 keV 
and added a Gaussian line near 3.2 keV to model a calibra- 
tion glitch l|Kubota et alJlsToiOh . We also included a narrow 
Gaussian line near 0.82 keV which was allowed to have sepa- 
rate normalisations between the front and back-illuminated 
units. Following Kubota et al., a 1 per cent systematic uncer- 
tainty was added to all XIS energy channels. PIN data were 
analysed from 18-50 keV, while using th e 'tuned' instru- 
mental background (|Fukazawa et al.ll2009h . The cosmic X- 
ray background (CXB) contribution to the PIN spectrunjf] is 
included as a spectral model componenlQ with 10 per cent 
freedom in its normalisation. We fitted for the normalisa- 
tion of the instrumental background spectrum, which has a 
nominal uncertainty of 3 per cen10. 

During the Suzaku observation, we obtained eleven 
RXTE pointings^ (ran ging from 1-11 ks apiece), which 
are interspersed throughout our Suzaku observing window 
and bracket it. We exclusively used the best-calibrated, 
'standard 2' spectra from PCU-2. These data improved the 
constraints on the continuum features, including the high- 
energy power-law component and the Compton hump. The 
RXTE spectra were background subtracted, corrected for 
detector dead time, and anal ysed from 2.55-45 keV. 

As in previous work (see lSteiner et aill2010l . [201ll ). we 
standardised the shape and normalisation of each detector's 
calibration so that the power-l a w par ameters of Crab spec- 
tra match the iToor fc Sewardl (|l974T l values; we introduce 
a floating cross-normalisation between detectors to account 
for any residual difference (< 5 per cent), except for the 
PIN which is assigned a fixe d normalisation o f 1.16 rel- 
ative to XIS - 0. M otivated by iTsuiimoto et all |201ll ) and 
llshida et al.1 (|201ll ). which show the back-illuminated XIS 
detector yields a significantly different spectral index than 
the two front-illuminated detectors, we include an extra pa- 
rameter for the difference between the XIS-1 spectral in- 
dex and those of XIS-0/3. This is incorporated into the 
model as a fit parameter and found to be quite modest: 
Ar ?s 0.015 ± 0.005 (Table [TJ. The additional freedom in 
the model significantly improves our fits (Ax 2 ~ 25), and 



4 http:/ /space. mit.edu/CXC/software/suzaku/acatt. html 

5 http:/ /space. mit.edu/CXC/software/suzaku/pest. html 

6 The CXB contributes negligibly to the XIS. 

7 http:/ /heasarc. nasa.gov/docs/suzaku/analysis/pin_cxb. html 

8 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

9 To avoid contam ination from the nearby PSR B0540-69 
l|Haardt et alj|200ll l. RXTE was offset by 0?25 in the opposite 
direction, and we have generated an off-axis response calibration 
to correct for this, while including a 1 per cent systematic uncer- 
tainty in each spectrum. 
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insignificantly alters the value of spin (the change is < 10 
per cent of the statistical uncertainty). 



3 MODELS 

We introduce two similar models, which we refer to below 
as Model 1 and Model 2, that differ in our treatment of the 
reflection component. Both models assume that this compo- 
nent is generated in the same manner: by disc photons that 
are Compton scattered in the corona and then reprocessed in 
the accretion disc. In t he first model, we employ the reflion 
|Ross fc Fabianl 120051 ) family of refl ection spectral models, 
and i n the second, use the pexrav (|Magdziarz fc Zdziarskil 
1 19951 ) family. In both cases, we mod e l the thermal disc emis- 
sion using KERRBB2 dLi et al.ll2005l ; iDavis fc Hubenvlliooil ; 
iMcClintock et ail 120061) and t he C ompton component via 



2009b. 2011). 



Simpl-R (jSteiner et al 

For the REFLiON-based model , we use refbhb 
(|Ross fc Fabianll2007l ; lReis et all 120081 ). Unlike its counter- 
part, reflionX, which describes reflection from supermas- 
sive BHs in active galactic nuclei (AGN), refbhb accounts 
for the Compton broadening and other effects unique to the 
hot, dense discs around stellar-mass BHs. The model ref- 
bhb includes a pre-packaged thermal component. This disc 
component is a single-temperature blackbody, and hence is 
intrinsically narrower than the standard multi-temperature 
disc spectrum. We remove this hardwired blackbody com- 
ponent by the following procedure: We construct a lookup 
table of blackbody parameters that match the refbhb in- 
puts, and we use this table to subtract off the best-fitting 
2-20 keV blackbody spectrum, while pairing the reflection 
portion of refbhb with the disc model kerrbb2. Although 
this approach does not account for radial variations in the 
reflection spectrum, it has the virtue of allowing refbhb to 
be applied more flexibly to spectra with a strong disc com- 
ponent. 

We henceforth refer to this custom version of refbhb as 
refbhb-M. The composite of refbhb-M operating jointly 
with Simpl-R and kerrbb2 comprises what we refer to here- 
after as Model 1. We also consider an alternate 'prime' ver- 
sion of Model 1 - Model lp - with refbhb-M replaced by its 
AGN counterpart reflionX. Apart from this substitution, 
Model 1 and Model lp are identical. 

For the reflection components of each of the models con- 
sidered here, we assume Solar metallicity, which is the only 
setting available for refbhb and refbhb-M. The LMC as a 
whole is known to have a lower metallicity; however, LMC 
X-l is likely t o be relatively me tal rich because of its young 
age (~5 Mvn lOrosz et alj|2009h . 

Our second model family uses ireflect which self- 
consistently computes the reflection edges for ionised gas, 
given an arbitrary input spectrum, ireflect is a generali- 
sation based on the model pexriv, but whereas pexriv is re- 
stricted to a pure power law for coronal emission, ireflect 
is freely combined with any coronal spectrum desired. As 
input to ireflect, we use the Compton component gener- 
ated in the accretion-disc corona, which (as above) is gener- 
ated by Simpl-R acting on the disc component kerrbb2. We 
hereafter refer to this composite model as Model 2. Although 
ireflect offers the advantage of computing reflection self- 
consistently given an arbitrary input spectrum for the illu- 



minating hard X-rays, it has several major drawbacks com- 
pared to refbhb/refbhb-M. For example, ireflect only 
computes edge absorption; that is, it neglects all fluores- 
cent line emission, including the Fe line itself. We insert the 
ionised Fe Ka line emission by adding an intrinsically nar- 
row Gaussian. Like refbhb-M, ireflect does not account 
for radial variation in the disc's structure. 

In summary, we face a tradeoff: Model 1 provides an 
optimal description of the atomic features, both lines and 
edges, whereas Model 2 provides a superior description of the 
'continuum' shape of the reflection component, especially at 
low energies. More specifically, although these two models 
describe the same process, Model 2 consistently describes 
the continuum shape determined jointly by the thermal, 
Compton, and reflection components, while fluorescent line 
emission is not incorporated self consistently. Conversely, in 
Models 1 and lp the atomic physics is self consistent, i.e., 
the treatment of emission and absorption is unified, but the 
shape and intensity of the reflection component is not di- 
rectly tied to the flux or curvature in the associated power- 
law component. Lacking an ultimate model of reflection in 
which the virtues of both models are captured, we employ 
both approaches separately in arriving at our final result. 

For all formulations (Model 1, Model lp, and Model 2), 
the inner-disc reflection is convolved with the relativistic 
smea ring kernel relconvf (|Dauser et al.ll201ol ; lFabian et al.l 
120121 ). We also include an unblurred reflection component 
to account for reflection far from the BH, which we demon- 
strate is produced by fluorescence of the stellar wind by the 
X-ray source (Section 15 . 3p . We model this sharp reflection 
component using reflionX. 

Photoelectric a bsorption is treated using tbvarabs 
llWilms et al. I I2000PI Our three model formulations in 
xspec notation are: 

( 1 ) tbvarabs x (simpl-r®kerrbb2+relconvf(g> 
refbhb-M+reflionX) , 

(lp) TBVARABS X (SIMPL-R(g>KERRBB2+RELCONVF® 

reflionX+reflionX) , 

(2) tbvarabs x (simpl-r®kerrbb2-|-relconvf® 
[ireflect(simpl-C)+gaussp c ] + reflionX) , 

where SIMPL-C is shorthand for the power-law emission iso- 
lated from SIMPL-R® KERRBB2 . 

3.1 Method 

Our results are obtained by first achievin g a set of pre - 
liminary spectral fits using xspec (v12.7; lArnaudl Il996l ). 
These fits provide the seed for a more intensive and ro- 
bust analysis using a Markov chain Monte Carlo (MCMC) 
ro utine, which is implemented via a package described 
in ISteiner fc McClintockl (120111 ). For this application, the 
MCMC routine has been modified to work with XSPEC. For 
each fit, a total of 3x10 s elements are generated"!, which 
is sufficient to reach convergencj^l. The results so obtained 



10 Th e ISM composition in TBVARABS is taken from lHanke et all 
j2009h . 

11 An additional fa 10 5 elements were generated during training 
and burn-in phases but these w ere not used in the final analysis 

sec ISteiner fc McClintocklboilJ). 

12 In all satisfactory iGelman fc Rubinl dl992h conver- 
gence diagnostic of R < 1.2 was obtained. 
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have been verified using xspec in conjunction with a se cond 
MCMC sampler: emcee l|Foreman-Mackev et alj|2012h . 

All but four free parameters were modelled with non- 
informative priors, either flat - for shape parameters such 
as r or T, or log-flat (i.e., a flat weighting on the log 
of the parameter) for scale parameters such as M and 
component normalisations. Two exceptions to this rule are 
Ah, which was assigned a normal prior di stributional of 
A(1.15, 0.15) x 10 22 cm" 2 (lHanke et al.ll2009T ), and the PIN 
instrumental background normalisation, which was taken 
to be a normal distribution of iV(l,0 .03) bounded withi n 
1 ± 0.06. Lastly, based on the results of lFabian et al.l (|20 12T ) . 
we impose a restriction on the paired va l ues of a* and the in- 
ner emissivity index qi. In lFabian et al.l (|2012T ). it was shown 
that for all non-maximal values of spin (a» < 0.94), qi is not 
expected to deviate substantially from its nominal thin-disc 
value of q — 3. High or extreme values of q\ are expected 
only for strong distor tions of spacetime near the horizon, 
for R < 2GM/c 2 (e.g.. IWilkins fc Fabianlboill ). To prevent 
runaway to unreasonable regions of parameter space, we as- 
sign a weakly constraining prior on q\ , a normal distribution 
of the following form: 

p{qi\a„) = N{3, 0.475 + 0.05 x (1 - a») _1 ). (1) 

Under this prior, at the minimum possible spin, a» = — 1, 
the prior's standard deviation is 0.5, and it increases slowly 
with spin, reaching unity at a* ~ 0.9. Thus, for most allowed 
values of spin, the fit is loosely constrained to conform to 
its nominal value, qi — 3. At the highest values of spin, the 
prior offers essentially no constraint because the standard 
deviation diverges. 



4 ANALYSIS 

Before applying these spectral models in detail, we first step 
back and adopt a more data-centric approach that firmly 
establishes the existence of a relativistically broadened Fe 
line. LMC X-l has a sibling BH binary in the Large Magel- 
lanic Cloud, LMC X-3. LMC X-3 is also persistently bright, 
and it is nearly always in a soft state. However, in contrast 
to LMC X-l, LMC X-3 is highly variable w ith L x /L E dd 
ranging from < 1 per ce nt to ~ 50 per cent l|Wilms et al.l 
l200ll : ISteiner et alj|20ich . Moreover, the spectrum of LMC 
X-3 app ears featureles s, to su ch extent that it was recently 
used by iKubota et alj |2010) as a benchmark for testing 
the performance of spectral disc models, which makes it an 
ideal reference source for our purposes. Use the ~ 70 ks 
Suzaku spectrum of LMC X-3 obtained in December 2008, 
studied by Kubota et al., and following precisely the same 
reduction procedure we used for LMC X-l, we compute the 
ratio of the spectra of the two sources for the combined 
front-illuminated XIS units. The virtue of this approach is 
that this ratio spectrum is completely free of any detector- 
calibration or data-reduction artefacts. 

Because the two continua have different shapes, the ra- 
tio spectrum has some smooth curvature, which we remove 
by fitting it to a low-order polynomial. The residuals clearly 

13 N([i, a) refers to a normal distribution centred upon fi with a 
variance of a 2 . 
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Figure 1. Top: Ratio of the spectrum of LMC X— 1 to that of 
LMC X— 3 shown fitted to a polynomial (solid line). Bottom: The 
residual spectrum, which clearly reveals the presence of broadened 
Fe Kct emission. 

reveal the presence of a broad Fe-line component in LMC 
X-l, as shown in Fig.[T] The line is peaked near 6.7 keV, and 
it has a broad red wing, as expected for a gravitationally- 
redshifted line fluoresced in the inner disc. 

The broad compone nt extends to bel ow 5 keV at an 
inclination of i = 36° (jOrosz et al.l |2009j) , which implies 
that the disc inner radius is r < 2.5GM/c 2 and the spin 
is correspondingly high. While one cannot obtain a defini- 
tive measurement of spin using this ratio spectrum alone, it 
nevertheless serves to confirm the presence of LMC X-l's 
broad Fe line and suggests that the BH is rapidly spinning. 

As a first fit to the spectral data, without yet treating 
effect of relativistic reflection, we apply a model of Comp- 
tonised disc emission combined with a reflection compo- 
nent which is unblurred and narrow (tbvarabs x (simpl- 
R®kerrbb2+reflionX)). The data-to-model ratio for the 
best fit (x 2 / ly 1-2) is shown with Fig. [21 in which the resid- 
ual features from 4-6 keV highlight an excess attributed 
to an unmodelled broad red wing of the Fe line. Likewise, 
the deficit above 8 keV and the excess and slight curvature 
above 20 keV indicate, respectively, the smeared Fe edge 
and Compton hump (although the Compton hump signal 
is statistically marginal). These residual features are indica- 
tors of an untreated component of relativistically broadened 
reflection. 

Having established the existence of gravitationally- 
broadened spectral reflection features, we will now apply 
the suite of disc/reflection spectral models described above 
and thereby constrain LMC X-l's spin. 

4.1 Spectral Analysis 

We begin by making the standard assumption of alignment 
between the BH and the binary orbital plane (i « 36°; 
lOrosz et al.ll2009r ). and we proceed to apply Models 1, lp, 
and 2 to the Suzaku data alone. We then follow by applying 
these same models to the combined Suzaku and RXTE data 
set. The dominance of the thermal component in the spec- 
trum of LMC X-l, and the modest luminosity of the source 
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Figure 2. Residuals to the best fit of a model with Comptonised 
disc emission and unblurred (narrow) reflection reveal a promi- 
nent broad red wing to the Fe line in the ratio spectrum. The PIN 
data show a hint of curvature in the energy range of the Compton 
hump. These features further indicate the presence of relativistic 
reflection from the inner disc. Here, the data have been rebinned 
for plotting purposes only. To allow for a straightforward com- 
parison, the inset highlights the same energy range as Fig. [T] 



(Section [I}, assure that the spectrum is well described by 
standard thin-disc theory. Therefore, in modelling the inner 
reflection component, we begin by assuming that illumina- 
tion of the disc follows a thin-disc p rofile with the stand ard 
value of the emissivity index, q = 3 (jFabian et al.lll989l ). 

Fig. [3] shows our best fits to the Suzaku spectra for 
all three models. The fitting results are summarised in Ta- 
ble Q] (q — 3 results are denoted by subtypes (i) and (hi)). 
MCMC is designed to directly provide posterior probability 
densities for the free parameters of our models. However, 
for large numbers of free parameters MCMC alone does 
not provide an optimised estimate of the minimum value 
of x 2 l v - Therefore, after running the MCMC chains, in or- 
der to provide an estimate of the usual goodness of fit, we 
have optimised the model \ 2 /v about the central values in 
the chain; these optimised fits yield the \ 2 l v values given in 
Table[U Although the MCMC chains do not directly deliver 
optimised goodness-of-fit estimates, we emphasise that the 
chains themselves provide the most direct estimates of the 
probability distributions for our model parameters, which is 
of chief interest. 

We explore the effect of allowing q to vary and to take on 
the form of a broken power law (model subtypes (ii) and (iv) 
in Table [1]). For moderate or low spins, such as those given 
for Model lp, this exercise is not physically motivated, but 
we nevertheless use this approach to understand the effect 
of having initially fixed q at its canonical value of 3. Once 
qi is freed, all three models return high values, qi > 3; the 
outer index is fixed to q2 — 3. As a net result of this exercise, 
the values of spin are slightly depressed. The dependence of 
a, on iji is illustrated in Fig. 2] There is generally a weak 
anti-correlation between the two parameters, which is most 
pronounced for Models lp (center) and 2 (right). 

Because a high value of the inner emissivity index qi in- 
dicates substantial light bending in the strong-gravity envi- 
ronment close to the black hole's event horizon, one expects 



high q\ to be associated with high values of spin. That we 
find the opposite correlation in our fits - high q\ is associated 
with low spin - is simply a systematic artefact of either our 
model or the data, which is possibly related to the limited 
signal-to-noise in the reflection features. 

As is evident in Fig. [3] the spin is obtained by decom- 
posing the blurred and sharp reflection components, with a 
substantial fraction of this signal coming from around the Fe 
line and edge. The sharp reflection is produced by material 
distant from the source, which we identify with the wind of 
the companion star. This subject is discussed in Section[5]3] 

Overall, the three models perform comparably well. 
(The statistical differences in the quality of fit between Mod- 
els 1, lp, and 2 are even less when one factors in the influ- 
ence of the priors - an effect not captured by \ 2 l v -) How- 
ever, because Model lp uses a reflection model designed for 
AGN, whereas Model 1 is optimised for the densities and 
environments of BH binaries, we favour Model 1 because it 
is physically more appropriate. Model 2, which offers the 
most self-consistent treatment of the broad reflection con- 
tinuum (while ignoring the fluorescent emission features), 
also returns a high spin. Model lp is less constraining, but 
consistent with the others. Although we adopt the results 
for Model 1 which we consider to be the most physical of 
the three models, we use the combined results for all the 
models in all the columns of Table [1] to account for sys- 
tematic uncertainty in our spin measurement. In particular, 
although their paired fits for qi and a* are physically sus- 
pect (described above), we nevertheless include Model lp(ii), 
Model 2 (ii) , and Model lp(iv) in determining the net result. 
This conservative approach increases our final uncertainty. 

Our final result for Model 1, combining all four variants 
(i.e., (i)-(iv) in Tabled]) is a* = 0.97±£;°3. We account for 
systematic effects by combining the results for all the models 
in all the columns of Table [T] (with equal weight assigned to 
each), and thereby obtain a„ = 0.97H*Ho'25. As a firm lower 
bound, at 99 per cent confidence, the spin is constrained 
a, > 0.2. This distribution, computed from the sum of all 
the MCMC chains, is shown in in Fig. [S] 

In order to assess the effect of assuming Solar abun- 
dances, we explore a reduced metallicity of Z — O.SZq 
for Model lp (this metallicity setting is possible using re- 
flionX but not refbhb-M). We find an increase of < 10 
per cent in the ISCO radius, which implies a reduction in 
the spin of Ao* w 0.02 from our fiducial values. Metallic- 
ity is, therefore, a relatively minor source of systematic un- 
certainty. Similarly, we have explored the uncertainty in- 
troduced by the 1.9 deg error in the orbital inclination 
l|Orosz et al.l 1200*91 ) and find that the resulting spin uncer- 
tainty is quite small (Aa* < 0.01) compared to our net un- 
certainty. Lastly, one notable deficiency of Model 2 is that it 
omits treatment of Compton broadening of the Fe line due 
to scattering in the (keV-temperature) disc atmosphere. Be- 
cause this temperature is far below the energy of the Fe line, 
Compton broadening produces an intrinsic red tail t o the 
line (e.g.. iTorreion et aljlioiol ; IWatanabe et al.ll2003l ). We 
explore this broadening by replacing the narrow Gaussian of 
Model 2 with a one-sided Gaussian with width a = 0.5 keV. 
The width used is approximately twice as large as that due 
to single scattering in the disc. The broadened line returns 
significantly improved fits (A% 2 > 15), which we interpret 
as a likely detection of some degree of Compton broaden- 
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Figure 3. The unfolded spectral fit and fit residuals for LMC X— 
1 using all three spectral models with q = 3. We show the fit for 
Model 1 (refbhb-M), our adopted model, in panel a; Model lp 
(reflionX) in panel b; and Model 2 (ireflect) in panel c. The re- 
flection signal is comprised of two parts: The primary comes from 
the inner portion of the BH's accretion disc (green lines), and 
the secondary from the distant gas in the system (purple lines). 
Comptonised disc emission is shown in red, and the composite 
model is overlaid in orange. In panel a, because REFBHB-M op- 
erates as reflection/reprocessing modifying a thermal signal, the 
relativistic reflection is drawn showing REFBHB-M+1.5/sc( SIMp L- 
R(giKERRBB2), arbitrarily chosen in order to depict a compara- 
ble normalisation to panels b and c. The XIS spectra have been 
binned to a minimum signal-to-noise of 70 for plotting only. 

ing. In all instances of Model 2 (i.e., (i)-(iv) in TableQ]), the 
spin returned is in the range a, — 0.94 — 0.98. We conclude 
that, because Compton broadening is an order of magni- 
tude smaller than the maximum gravitational redshift for 
a rapidly-spinning source, the impact on the fitted spin is 
moderate. Generally, accounting for this extra broadening 
causes the model to return a slightly lower value for spin in 
Model 2. 



5 DISCUSSION 

5.1 The Connection Between the Compton Power 
Law and Reflection 

We now explore the relationship between the most promi- 
nent feature of the reflection component, namely the Fe 
line, and the degree of Comptonisation. To do this em- 
pirically, we make use o f a si mplified continuum model 
employed in iGou et alJ l|2009l ). with a small modifica- 
tion: TBVARABS(SIMPL®KER RBB2+GAU SSFfif'^l. We apply 

this model to both the full IGou et all data sample of 53 
RXTE spectra and to the RXTE spectra collected during 
our Suzaku observation. This model lacks the physical rigour 
of those we have been considering, but it has the virtue of 

14 For all spectra, we use fixed values for the Gaussian line energy 
and width (6.5 keV and 0.25 keV respectively), and we use the 
continuum-fitting value of spin obtained by Gou et al. (a* = 
0.92). Further following IGou et al.l (2009), we restrict the power- 
law index to values between T = 2.25 and V = 3.25. 



Figure 4. Selected correlation results from the MCMC analysis 
showing that spin and the inner emissivity index qi are weakly 
correlated. Red contours show lcr, 2a, and 3<r confidence limits 
about maximum density values. Top panels correspond to Suzaku 
fitted alone (model subtype (ii)), and the bottom panels show re- 
sults for RXTE and Suzaku fitted together (model subtype (iv)). 
Panels on the far left use REFBHB-M (Model 1), those in the cen- 
tre column use reflionX (Model lp), and those on the right use 
IREFLECT (Model 2). 
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Figure 5. Our final and adopted spin probability distribution 
obtained by combining results from all the models in Table [l] 
with equal weights. 

being sensitive to the strengths of the Fe line and Compton 
component without requiring additional parameterisatiorF^l. 

As shown in Fig. [6l a clear and positive correlation is 
observed which is well-described by a linear fit (x 2 jv = 1.17) 
with a correlation that is significant at > 12cr. Regardless of 
the geometry of the corona - whether a disc-hugging skin, 
a lamp-post, or a centrally-concentrated cloud - the exis- 
tence of this correlation confirms that the coronal Comp- 

15 Because this simplistic model does not account for photoab- 
sorption of Compton-backscattered photons, which return to the 
disc, the value of /gc is lower here by a factor of roughly two 
compared to its value under Models 1, lp, and 2 above. 
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Figure 6. A strong correlation exists between the relative 
strength of the Compton power-law component (/sc) an d the 
Fe-line flux based on RXTE data. Because the simplistic model 
employed assumes an arbitrary line wid th, the line flux is shown 
to arbitrary scale. The iGou et al.l j 20091) spectra are displayed as 
closed symbols, and our observations are plotted as open symbols. 
A linear fit is overlaid as a dashed line. 



ton emission is the primary source that illuminates the disc 
and produces the reflection/fluorescence component in LMC 
X-l. This result indicates that disc self- irradiation (which 
could be important in producing reflection features either at 
high luminosities or in Compton-weak sources) is not domi- 
nant here. Our f inding s are in close agreement with those of 
Iffiemstra et all (|201lh . who demonstrate for the BH XTE 
J1652-453 the coexistence of a bright disc and strong re- 
flection features. Hiemstra et al. point out that strong coro- 
nal flows in bright (intermediate-state) discs can produce 
the most prominent broad Fe lines, which may explain the 
strong line feature we observe in LMC X-l. 



5.2 The Thermal Continuum 

It is often problematic to use the same spectrum to estimate 
spin by both the continuum-fitting and Fe-line methods be- 
cause the optimal spectrum for continuum-fitting has a weak 
power-law/reflection component, whereas a spectrum that 
is optimal for the reflection method generally has a domi- 
nant power-law component, which compromises continuum- 
fitting measurements. Unfortunately, for the Suzaku and 
RXTE spectra treated here, one cannot obtain a reliable spin 
measurement via continuum fitting. Despite its stable lumi- 
nosity, LMC X-l exhibits strong rms power fluctuation^^]- 
To obtain reliable continuum-fitting results, on e generally 
restr i cts the rms power to be < 8 p er cent (e.g.. IGou et al.l 
120091 ; iRemillard fe McClintockfooOfj ) . However, for the spec- 
trum in question, the rms power is in the range 27-35 per 
cent, and therefore a reliable continuum-fitting spin result 
cannot be obtained from these data. 



16 We follow the convention of IRemillard fe McCIintockl d2006f) , 
who define rms variability as the average power from 0.1-10 Hz 
in the band 2-20 keV. 



In a study of LM C X-l's variable power-law component, 
iRuhlen et all (|201ll ) examine the archive of RXTE observa- 
tions of LMC X-l, and, assuming a uniform and nominal 
power-law index, find a weak anticorrelation between disc 
flux and temperature. This is exactly analogous to the an- 
ticorre lation of spin and M shown in figure 4 of IGou et al.l 
(2009). Ruhlen et al. interpret their flux, temperature anti- 
correlation as implying the presence of a variable corona that 
sometimes obscures the inner disc. While this interpretation 
is possible, it implies the occasional presence of an optically- 
thick, ~ 50 keV coronal cloud, which has never been ob- 
served. We therefore consider this explanation unlikely. This 
anticorrelation may instead be related to instabilities in the 
disc, as suggested by the correlation between the strength 
of the intense power-law co mponent and th e rms variabil- 
ity in the power spectrum |Gou et al.l 120091 : and with the 
reflection strengt h, see Sectionl5.l| ). Speculatively, the vari- 
ability noted by IRuhlen et al. | may be consistent with the 
inhomogeneous-disc model of Dexter fe Quatacrt l|2012ri. 

For any of these interpretations, the Gou et all ( 20091 ) 
continuum-fitting spin measurement is robust to these ef- 
fects because for measuring spin, Gou et al. selected only 
those data in which the Compton component is weakest 
(/sc < 0.075), or equivalently, data for which the rms vari- 
ability is weak. However, the continuum-fitting measure- 
ment is confounded for the Suzaku observation considered 
here. 

We note that the sign and magnitude of the deviation 
we find for t he continuum-fitt ing value of spin compared to 
the result of IGou et all <|2009D is consistent with what has 
been obtained for other systems; namely at excessively large 
values of rms, the continuum-fitting model te nds to return 
errone ously small values of spin (e.g., fig. 1 in lSteiner et al.l 
l2009al ). 

Because in this instance the ther mal continuum mode l 
is flawed, and additionally because iKubota et all (|2010h 
emphasise that a true disc spectrum may be intrinsically 
broader than the zero-torque kerrbb model, we examine 
the effect of using a modified torque at the inner-boundary. 
For the reasonable values we have explored, 77 < 10 per cent 
(see iLi et all 120051 for details on the torque prescription), 
the reflection spin measurement is insignificantly affected 
by this change in the continuum model. To otherwise assess 
the impact of our choice of thermal contin uum model, we 
have tested replacing kerrbb2 with diskbb jMitsuda et all 
1 19841 ). and also with bhspec (|Davis fe Hubenvl 120061 ) in 
Models l(i), lp(i), and 2(i). The goodness-of-fits returned 
with these disc models is worse by A\ 2 ~ 20, and the spin 
in each instance is consistent with a, > 0.9. We conclude 
that our reflection spin result is robust to the choice of ther- 
mal continuum model. 



5.3 Distant Reflection 

Based on optical in tegral-field s p ectros copy with a res- 
olution of ~ 1 pc, ICooke et all (|2008l ) mapped a large 
(~ 4 pc) cone-shaped ionisation nebula enshrouding LMC 
X-l (|Pakull fe Angebaultl Il98rj ). The cone has an open- 
ing angle ~ 45° and a mean density on large scales of 
n < 100 cm -3 . Meanwhile, our spectral fits imply the pres- 
ence of a strong source of reflection that is distant from 



8 Steiner et al. 



100 r 



80 - 



60 - 



S 40 



20 - 




1.5 



- 1.0 



0.5 



J 0.0 



12 3 4 

Radius (pc) 

Figure 7. Radial profiles for the density and gas temperature 
in the ionisation cone are approximately flat. This supports our 
assertion that the optical ionisation cone is a shell of shocked 
interstellar gas. This gas encloses a bubble inflated by the stellar 
wind (e.g., compare with the 'shell' region in fig. 3 of lWeaver et al.l 
ll977lV 



LMC X-l (i.e., separated by r >> 10 cm) with a luminos- 
ity ~ 10 36 - 10 37 erg s" 1 . 

Both the cone structure and the reflection component 
are readily explained at once by positing the presence of a 
powerful wind that is photoionised by the luminous X-ray 
source. We show that the wind is likely supplie d by a domi- 
nant outflow from the O supergiant companion (jOrosz et al.l 
2009) (possibly coupled with a fast wind launched from the 
disc). The source of illumination, meanwhile, is the self- 
irradiated accretion disc, which supplies the requisite flux 
of ionising photons after reprocessing < 5 per cent of the X- 
rays from the inner disc i nto ionising photons in the outer 
disc (radii > WOGM/c 2 ; [Gierliriski et "ail 120081 : IcTou et al.1 
I2009T ). 

Adopting a very simplified model in which reflection is 
approximated as efficient Thomson scattering through an 
idealised wind with ionisation parameter ~ 10 3 erg cm s — 1 , 
we find that the wind is launched at a distance ro ~ 0.2 AU 
from the BH and with a density of no ~ 10 10 cm" 3 . These 
values match the BH-to-compan ion star distance and the 
surface density in its stellar wind (jLamers fc Leithererll 19931 ; 
lOrosz et al.ll2009l ); this identifies the stellar wind as the 
likely source of distant ionised gas. At sub-pc scales, the 
ionisation of the wind is constant (n oc r~ 2 ), but because 
density falls off rapidly, most of the (sharp) reflection is cen- 
trally concentrated and emitted from the innermost several 
AU. 

Meanwhile, at parsec scales, the wind is shocked against 
the ambient interstellar medium (ISM). As a result, the ioni- 
sation cone is mostly filled with a dense gas of shocked wind. 
This cone, in turn , is enclosed in a th in shell of even denser 
shocked ISM (e.g. JWeaver et aT]|l977l ). This picture is borne 



out by the flat temperature and density profiles!] within the 
ionisation cone, which are shown in Fig. [7] 

Adopting this simplistic description of the wind as valid, 
and employing the I Weaver et all l|l977f ) wind-bubble model, 
the projected mass in the wind is of order ~ 0.1M©, which 
suggests a time-scale of ~ 10 6 years to produce the ionisa- 
tion cone surrounding LMC X— 10 The mass in the cone is 
then commensurate with the integrated mass loss from the 
companion star, and a time-scale resu lts which is compa ra- 
ble to the » 5 Myr age of the system l|Orosz et alj|2009h . 



6 CONCLUSIONS 

We have analysed a deep Suzaku and RXTE observation of 
LMC X-l, and we have demonstrated that a broad Fe Ka 
line is present in the spectrum of the source. Using existing 
and new spectral reflection models, we have measured the 
BH's spin parameter to be a* = 0.971q'^3 for our favoured 
model (Model 1), and to be a* = 0.971q 25 when making an 
allowance for systematic error by considering our full range 
of models. At 99% confidence, we establish a, > 0.2. Both 
of these spin estimates are in agreement with the spin deter- 
mine d using the X-ray continuum-fitting method (|Gou et al.l 
l2009h . 

Apart from the measurement of spin, we present two 
additional results: (1) For a large sample of RXTE spectra, 
we demonstrate a strong and positive correlation between 
the Compton and reflection components. (2) Far from the 
strong gravity environment of the BH, we identify the wind 
of the massive companion star as the source of a luminous 
and sharp reflection component. We conclude that the wind 
and persistent X-ray source together maintain the parsec- 
scale ionisation cone that envelops the binary system. 
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17 These profiles are derived from the data in lCooke et al.l d2008h 
by using the [S II] AA6716/6731 ratio as a density diagnostic and 
the [O III] AA(4959+5007)/43 63 rat io as a tempe r ature diagnostic 
dOsterbrock fe Ferlandl 12006^ . See ICooke et al.l j200Sh for more 
details. 

18 Derived using the ISM density in LMC X— l's local environ- 
ment (an O-B association) of ~ 20 cm -3 . This, in turn, is ob- 
tained by assuming the standard jump condition at the outer 
shock boundary given that the density in Fig. [7] corresponds to 
the high-density outer shell. 
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!000 +iooo 

' uuu -2000 



960+ "" 1 
U.»DU_ 003 

0.919 ± 0.003 

-0.014 ± 0.003 

8.9 ± 0.7 

1.04 ± 0.01 

1.16 



1 41+0-01 
-0.03 

2.62 ± 0.05 
0.13 ± 0.01 
0.82 ± 0.01 

1.63 ± 0.02 
4.5+ 



S + 0.4 
-0.2 



4.9^ 



Q + 0-1 

0.82+|^ 



14000+ 2000 



960+ n 001 
U.»DU_ 003 

0.918 ± 0.003 

-0 015+ ' 002 
u.uio_ 003 

8.9 ± 0.7 

1.03 ± 0.01 

1.16 



x 2 /" 


395.4/396 


396.3/398 


415.0/395 


389.0/394 


375.3/396 


387.3/393 


Parameters 


Model l(iii) 


Model lp(iii) 


Model 2(iii) 


Model l(iv) 


Model lp(iv) 


Model 2(iv) 


Instruments 


Suzaku 


Suzaku 


Suzaku 


Suzaku 


Suzaku 


Suzaku 




& RXTE 


k RXTE 


& RXTE 


& RXTE 


& RXTE 


& RXTE 



JV /f (10 22 cm -2 ) 

r 

/sc a 

«.,CF b 
M{l0 18 g s" 1 ) 
91 ° 

«br(«ISCO) C 



-0.005 



1.3643 

2.66 ± 0.03 

0.183 ± 0.007 

n 8K1+ 006 
0.8bl_o 004 

1.79 ± 0.02 
3 



3 ) 



-Hdcn (10 cm" 
kT Iof (kcV) 
Illum/BB 

JV rof (10" 7 )/(10- 3 ) d 
£ c (crg cm s ) 
T rcf (10 6 K) 
E Fe {keV) 

F Fc fl0~ 4 ph s^cm" 2 ) 
?D°(8rg cm s" 1 ) 
iV rcf , D flO- 7 ) 
^xisa/^xiso 
'Vxisi /-^xiso 1 
Ar(XISl - XIS0)' 

N PIK CXBiiO" 4 ?!! s~ 1 cm~ 
PINBG.Inst Norm k 

iVpi N /iVxiso h 
^rxte/^xiso 11 



0.971 
6.4 4 



7 +0.01 
' -0.05 
4 +0.4 
*— 1.2 

u - M <-0.007 

034+""" 1 
U.UJ-±_ 002 

3.3 ± 0.3 



5000_ 1500 . 
50+ - 14 
0.960 ± 0.002 
915+ 002 
-0.018 ± 0.001 
9 4+ 02 
1.03 ± 0.01 
1.16 
0.999 ± 0.011 



1.55 ± 0.02 
2 74 +0.04 
2 ' '*-0.05 
0.21 ± 0.01 
„ +0.008 
U - 8ii -0.00G 
. 77 +0.03 
L - ' ' -0.03 
3 



, + 0.06 
-0.27 



1.1 ± 0.2 

9500+ 200 



600 ± 200 
8 ± 4 
0.960 ± 0.002 
0.918 ± 0.004 
-0.015 ± 0.003 
9.1 ± 0.6 
1-06-0.02 
1.16 
989+ 006 



1.44 ± 0.01 
2 7n +0.01 
2 -' u -0.03 

143+°"" :! 

0.810 ± 0.005 
1.65 ± 0.02 
3 



OQ2 + -003 
-0.060 



7400+ 1100 
8.9 ± 0.2 



1.39 ± 0.01 

2.66 ± 0.03 

n 1S n+ 009 
o.iau_ oo 6 

0.857 ± 0.005 

1.81 ± 0.03 

4 4+ 02 
-0.4 

3.0 ± 0.5 

n S4+ - 03 
0-84_ 07 

5.3 ± 0.3 

0.741 ± 0.007 

0.032 ± 0.002 

3 4+ 02 
J ^-0.3 



1.58 ± 0.02 
2 7 g+0.04 
2 - ,D -0.02 
0.210 ± 0.007 
0.805 ± 0.006 
i-»u_o 02 
4.4+" " 



,+0.3 
-0.4 
4.3 ± 0.5 
, +0.10 
-0.14 



0.41] 



6.05 



+ 0.03 



1.0 T 



n+0.1 
-0.4 
900 ± 200 
,+ 19 



-0.6 



3.3: 

0.959 ± 0.002 

0.917 ± 0.001 

-0.014 ± 0.001 

9 5+ 03 

1 041+ 004 
r.u I ±i_ oo6 

1.16 

1 049+ 005 
1^^^-0.004 



5000_ 1000 
0.5 ± 0.1 
0.961 ± 0.002 
0.918_ ooi 
-0.015 ± 0.002 
9.4 ± 0.2 

1 041+ 004 
-0.009 

1.16 
1.001 ± 0.011 



i q+O-S 
5800 ± 1000 



540+ 50 
° u -150. 

12 ± 1 

0.960 ± 0.002 

916+ 005 
u.yro_ 003 

-0.017 ± 0.003 
« 1+0-5 

1.05 ± 0.01 

1.16 

n no-,+0.006 

0.983_Q 008 



1.48 ± 0.02 
2.76 ± 0.03 

0-108_ 007 

0.791 ± 0.007 

1.69 ± 0.02 

4 9+ 02 
*' a -0.3 

3.7 ± 0.4 

94+ - 01 
u - y ^_0.08 



4100 ± 1500 

6 66 +0 - io 
°- DD -o.oi 

1.3+" 



12011 



+ 0.4 
0.3 
+ 200 
100 



4±! 
0.960 ± 0.002 

918+ 001 
u.3±o_Q 002 

-0.016 ± 0.001 

9 8+ 02 

1 039+"'"" 6 
i.uo»_ 009 

1.16 

1.060 ± 0.005 



898.6/1121 



887.2/1123 



984.6/1131 



885.5/1119 



862.1/1121 



944.4/1129 



Note. — MCMC fit results. Uncertainties arc the minimum- width 68 per cent confidence intervals about the posterior mode from the 
MCMC runs (computed in the log for scale parameters), and X 2 / 17 values are listed for the best fit obtained. Numbers with no uncertainty 
have been set to fixed values. We note that there are unshown fit parameters for each of the RXTE spectral fits (Models l(iii) through 
2(iv)). In each of those cases, the eleven RXTE spectra arc each fitted for an independent value of /sc- Additionally, for Models l(iii), 
lp(iii), and l(iv) and lp(iv), each RXTE spectrum is fitted for its own value of iV rc f . These values arc distributed approximately as 
/SC = 0.19 ± 0.06 and JV rof = (3 ± 1.5) X 10 -3 for Model l(iii) and (iv); / sc = 0.17 ± 0.05 and JV rcf = (2 ± 2) X 10~ 7 for Model lp(iii) 
and (iv); / S c = 0.15 ± 0.06 for Model 2(iii) and (iv). 

a Thc fraction of disc photons scattered into the Compton power law. 

b Continuum- fitting parameters: the spin and mass accretion rate. Mass, inclination, and distance have been frozen at their nominal 
values from lGou et alj i2009l h See lOrosz et al.l i2009l ). 

c Rclativistic smearing indexes qi and q2 describe, respectively, the inner and outer power-law illumination pattern of the disc. The 
power law is broken between the two regimes at i?v, r ; 92 i s se f f° 3 ln an instances. 

d Normalisation for REFLIONX, in units of 10~ 7 , and for REFBHB-M, in units of 10 — 3 . 

c The ionisation parameter of the reflecting material. 

^Fe-line parameters i?p c and Fp B give the line energy and line flux, respectively. 

s The subscript 'D' refers to reflection emission distant from the BH, and hence not relativistically smeared (see Section l5.3l . 



